Integrable deformations of Lotka-Volterra systems 



Angel Ballesteros, Alfonso Blasco and Fabio Musso 

Departamento de Fisica, Universidad de Burgos, 09001 Burgos, Spain 
e-mail: angelb@ubu.es, ablasco@ubu.es, fmusso@ubu.es 



Abstract 

The Hamiltonian structure of a class of three-dimensional (3D) Lotka-Volterra (LV) equations is 
revisited from a novel point of view by showing that the quadratic Poisson structure underlying 
its integrability structure is just a real three-dimensional Poisson-Lie group. As a consequence, 
the Poisson coalgebra map A^ 2 ) that is given by the group multiplication provides the keystone 
for the explicit construction of a new family of 3N-dimensional integrable systems that, under 
certain constraints, contain N sets of deformed versions of the 3D LV equations. Moreover, 
by considering the most generic Poisson-Lie structure on this group, a new two-parametric 
integrable perturbation of the 3D LV system through polynomial and rational perturbation 
terms is explicitly found. 
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1 Introduction 



The Lot ka- Volt err a (LV) system of differential equations 



Xi — Xi 



,L (1.1) 



where atj and d{ are constants, were initially introduced for L = 2 in order to model prey-predator 
and chemical reaction dynamics [TJ [2] , but nowadays it is well known that the generic system (II. ID 
plays an outstanding role in many different nonlinear models with a wide range of applications (see, 
for instance [3] H] and references therein) . 

In particular, the three-dimensional (L = 3) LV system has attracted much attention, in many 
cases focused in the search for the specific values of the parameters aij and di for which a constant 
of the motion does exist (see [H EJ El [U EJ [10], [H] ) . Obviously, such integrability is guaranteed in 
case that the LV system could be obtained as the Hamilton equations x = {x,l-L} provided by a 
Hamiltonian H and a suitable three-dimensional (generalized) Poisson structure |12[ H3| [141 [T5| [16] 
[T71 118]. In particular, in [12] it was shown that an integrable LV system previously described in [7] 
and given by 

x = x(Cy + z + A) 

y = y(x + Az + fi) (1.2) 
i = z(Bx + y + v) 

where 

ABC = -1 and v = /j,B — XAB (1.3) 
was bi-Hamiltonian. One of the Hamiltonian functions leading to (jl.2p is given by 

% = ABx + y-Az + u In y — (J, In z (1.4) 

together with the quadratic Poisson structure 

{x,y} = Cxy {x,z} = BCxz {y,z} = —yz (1.5) 

which has the following Casimir function (that is indeed a constant of the motion for the LV 
system): 

C = AB\nx-B\uy + \nz. (1.6) 

Moreover, the second Hamiltonian structure for the same system is obtained by taking C as the 
Hamiltonian and by considering a different cubic Poisson bracket, for which T~i is just the Casimir 
function (see [12] for details). 

The aim of this paper is to generalize the abovementioned quadratic Hamiltonian structure of 
the LV equations by showing that it can be interpreted as a Poisson-Lie bracket on a certain Lie 
group. Therefore, these 3D LV equations can be understood as an instance of integrable Poisson-Lie 
dynamics, and this property will enable us to get new integrable deformations of the LV equations 
in two different ways. Firstly, by using the group multiplication (which defines by construction a 
Poisson map A^ 2 ) between two copies of the group manifold) in order to get an embedding of certain 
integrable perturbations of LV equations into a higher dimensional integrable system. Secondly, by 
looking for deformations of the underlying quadratic Poisson-Lie bracket that are preserved under 
the group multiplication. 

The structure of the letter is as follows. In the next section we introduce a novel Poisson- 
Lie group approach to the quadratic Hamiltonian structure of LV equations. In Section 3 we 
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use the group multiplication map A^ 2 ) in order to obtain new integrable perturbations of the 3D 
LV equations through the embedding method. In Section 4 the most generic deformation of the 
quadratic Poisson-Lie structure is found, and a new two-parametric integrable deformation of the 
LV equations containing polynomial and rational perturbation terms is explicitly obtained. A final 
section summarizing the main results and pointing out several open questions closes the paper. 



2 Lotka-Volterra equations as Poisson-Lie dynamics 

Let us consider the following quadratic Poisson algebra V |1U|, [TT] that generalizes (|1.5p 

{X,Y} = aXY {X,Z}=(3XZ {Y,Z} = ~/YZ. (2.1) 

The Casimir function for V reads 

C = X^Y p Z- a . (2.2) 
If we consider the Hamiltonian function -which also generalizes (|1.4p - 

H = ai X + a 2 Y + a 3 Z + b x \ogX + b 2 logY + 6 3 logZ (2.3) 
we get the following LV equations as the integrable dynamical system given by F = {F,7i}: 

X = X[aa 2 Y + (3a 3 Z+{ab 2 + f3b 3 )] 

Y = Yl-a^X + ja 3 Z + (763-061)] (2.4) 
Z = Z[-p ai X - 1C l 2 Y -{Pb l + 1 b 2 )\. 

These are Lotka-Volterra equations that depend on 9 free parameters, and the system is always 
integrable since the Hamiltonian and the Casimir function are constants of the motion. Note that 
if the three en parameters are all of them different from zero, then they can be reabsorbed in the 
Hamiltonian and equations through a linear change of variables, and the following skew-symmetric 
form of the LV system is obtained (see [TU] [TT]) 



X = X[aY + (3Z + (ab 2 + 0b 3 )} 

Y = Yi-aX + ^Z + ijbs-ah)} (2.5) 
Z = Z[-/3X - 7 Y -^h+ 7 b 2 )}. 

Moreover, this system can be also shown to be equivalent (through a further linear change of 
variables) to the more usual 6-parameter LV equations given by (jl.2|) (the so-called ABC-system). 
In fact, if a, f3, 7 ^ and we define 

111 , . 

ai = , a 2 = , a 3 = - (2.6) 

a 7 p 

then from (|2.4|) we get the ABC system with 

A-l b — — r--- 

A ~ P' B ~a C ~ 7 (2-7) 

\ = ab 2 + /3b 3 , = — a&i + 763, v = -/3bi — jb 2 

which fulfils the conditions 

ABC = -1, u = fiB-XAB (2.8) 
that characterize the bi-Hamiltonian case described in 1121. 



At this point it is important to stress that, as we shall see in Section 3, after deforming the 
Poisson bracket (I2.ip the abovementioned linear changes of variables that relate the different pre- 
sentations of LV systems will no longer work. Therefore in the rest of the paper we will consider 
as the LV system the 9-parameter fully generic case (12. 4p . 



3 



2.1 The LV Poisson algebra as a Poisson-Lie group 

Let us consider the Lie group G generated by the 3D multiparametric Lie algebra g a ,i3^ given by 



T 

[eo,eij = -ei 



[eo,e 2 J 



a 



e2 [ei,e 2 ] = 0. 



(2.9) 



A three-dimensional representation of this Lie algebra is given by: 



P(ei) 




p(e 2 ) 




p{eo) 



I 

o -I o 





(2.10) 



In this representation, the three-dimensional form of a generic Lie group element M reads: 



M = exp (E 1 p(e 1 )) exp (E 2 p(e 2 )) exp (E p(e )) 



exp(i J B ) 





E 1 

exp(-i^o) ^2 
1 



(2.11) 



By changing to the new group coordinates X = exp(Eo), Y = Ex, Z = E 2 , we obtain 





( 


X~? 


Y 


M = 




X~i 


Z 




\ 





1 



and the matrix product of two different group elements M 2 and M\ will be 



Mo-M 



xi 



Y 2 

_ 7 



x 2 
V o o 



z 2 
1 



xf 

x 1 
V o o 



o n 

_ 7 



l 



(2.12) 



V 










(X 2 X l ) 




X 2 P Y l + Y 2 

x 2 « z x + z 2 
1 



(2.13) 

By following the usual convention in quantum group theory (see |19j for details) , this group mul- 
tiplication can be rewritten in more algebraic terms by identifying the two sets of group coordinates 
{X\,Yi, Zi} and {X 2 ,Y 2 , Z 2 } with the 'tensor product variables' 



X 1 =X®1, 
Y 1 = Y <g> 1, 
Z x = Z ® 1, 



X 2 = 1 ® X 
Y 2 = 1 ® Y 
Z 2 = l®Z 



(2.14) 



that come from considering each set of group coordinates as a set of dynamical variables living on 
a different copy of the group manifold. Under this viewpoint, the product of two group elements 
reads 



M 2 ■ Mi 



( 



\ 



X~H ®Xt> 





Y ®X? + 1®Y 

X~i <g> X~i Z (g> X~i + \®Z 
1(g) 1 



(2.15) 



This algebraic approach suggests that the matrix elements of the product fj2. 15f) of two group 
elements can be interpreted as the images of a mapping A^ 2 ) : M — > M <S> M, namely 



[A( 2 )(x)]i o A( 2 )(y) 

M-2 • Mi - 1 [A( 2 )(X)]-i A( 2 )(Z) 

A( 2 )(l) 



(2.16) 
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Such mapping is called the coproduct (or comultiplication map) and reads 

A (2) (X) = X®X = X 2 X X 

A( 2 )(y) =Y®X? + 1®Y = XjY x + Y 2 (2.17) 

A (2) (Z) =Z<g>X~« + 1® Z = X~"Z 1 + Z 2 , 

besides A^ 2 )(l) = 1 ® 1. Obviously, this construction can be generalized to the product of N group 
elements, and then we get the explict form of the iV-th coproduct A^ N \ that maps M into the 
tensor product of N copies of M. 

Now, the key observation is the following invariance property: If we assume that the 'abstract' 
commutative algebra of smooth functions on the group coordinates {X, Y, Z} is endowed with the 
LV Poisson structure V given by (|2.ip . then it turns out that the coproduct A^ 2 ) is a Poisson 
algebra homomorphism with respect to V, i.e., the following relations hold: 

{ A< 2 > (X) , A< 2 > (Y)} = a A& (X) A< 2 > (Y) 

{A^(X),A {2 \z)} = p A^ {X) A {2 \Z) (2.18) 
{a( 2 )(y),a( 2 >(z)} = 7 A^(y)A( 2 )(Z), 

where in f|2. 18|) {, } denotes the natural Poisson structure on V ® V given by 

{a ® b, c ® d} = {a, c} ® bd + ac ® {b, d}. ( 2 -19) 

The previous statement is proven by a straightforward computation and means that (C ao (G),V) 
is a Poisson-Lie group (equivalently, that (V, A^ 2 )) is a Poisson coalgebra [191 120] ) . Therefore, from 
this perspective we can conclude that the LV equations (I2.4p are just a specific example of Poisson- 
Lie dynamics defined on (C°° (G) , V) and generated by the Hamiltonian 

H = a 1 X + a 2 Y + a 3 Z + b 1 logX + b 2 logY + b 3 log 2. (2.20) 

Therefore, this result endows the known quadratic LV bracket (12. ID with a new underlying Poisson- 
Lie group structure, and the coproduct map A^ 2 ) will be essential in order to get integrable defor- 
mations of the LV system. 



3 Integrable deformations induced from the coproduct map 

In the sequel we show that the Poisson-Lie nature of the bracket V provides a straightforward 
deformation approach to LV equations which is based on their embedding within higher dimensional 
integrable systems. These systems will be obtained by taking into account that Poisson algebras 
endowed with a coproduct map A^ 2 ) give rise to a systematic way of constructing integrable systems 
with an arbitrary number of degrees of freedom (for a detailed exposition of such coalgebra approach 
to integrability, including its relation with other approaches to integrability and many explicit 
examples, see [201 Elj l22| 123] ) . 

In order to make this construction explicit, let us work out the N = 2 case by considering as 
the new Hamiltonian the coproduct A^ 2 ) of the N = 1 Hamiltonian (12.3p . namely 

= oi A^ (X) + a 2 A< 2 > (Y) + a 3 A® (Z) + h A< 2 > (log X) + b 2 A< 2 > (log Y) + b 3 A^ (log Z) 
= a x (X 2 X 1 ) + a 2 (x}y x + Y 2 ^j + a 3 [x~^Z x + Z 2 ) 

+h log {X 2 X X ) + b 2 log (Xf Yx + Y 2 \ + 6 3 log (X 2 « Z x + Z 2 ) . (3.1) 



5 



Now, by considering the natural definition (|2.19p of the Poisson bracket on two copies of the group 
manifold, we get the following Hamilton equations for this 6D dynamical system: 



X 1 = 

Yi = 

Zi = 

x 2 = 

Y = 



Xi (aa 2 Y l xj +0a 3 Z 1 X^") + X l 
Yi {-aaiX x X 2 + 7a 3 ZiX 2 " " ) + yYi 
Z t (-paiXiX 2 - 1 YiX^ -7Z1 
X 2 (aa 2 Y 2 + /3a 3 Z 2 ) + X 2 



aYi 



b 2 x£ 



L 2 21 -f 12 

baZiX'- 



pZi 



b 3 X 2 



X%Y 1 +Y 2 J \X 2 «Z 1 + Z 2/ 

-abiYi 



\_X 2 a Z\ + z 2 
b 2 YiXf 



aY 2 



X*Y 1 + Y 2 \ 



Y 2 (~aaiX 2 Xi + ^a 3 Z 2 ) + ^Y 2 



XiYi+Y 2 , 
b 3 Z 2 



-PhZi 



(3Z 2 



X 2 a Z\ + z 2 



Xn " Zi 



-aY 2 



h + %x}Yi a 2 



Z 2 (-(3a 1 X 2 X 1 - -ya 2 Y 2 ) - -yZ 2 



X?Yi+Y 2/ 



— —X 2 " Zi I a 3 H — — 

a V Xn " Zi + Z 2 



b 2 Y 2 



L X 2 Yi + Y 2 



-f3Z 2 



h + ^X^Yi ia 2 



X£Yi+Y 2 



—X 2 a Zi ( a 3 
a \ 



X 2 a Zi + Z 2 



(3.2) 



This system is completely integrable since by making use of the Poisson coalgebra approach it 
is immediate to show [2D] that, besides the Hamiltonian T~L^ 2 \ there exist three more integrals of 
the motion in involution. They are the two Casimir functions for each subset of group variables 
together with the coproduct of the Casimir function, namely: 

d = X^Yfz^ 



C (2) . = A (2) (C) = ( Xl X 2 )^ X$Yl + Y 2 ) (X^Zi+Z; 



C 2 = XpYfz^ 



(3.3) 
(3-4) 



It is straightforward to realize that the previous six equations are two 'twisted' sets of integrable 
deformations of the 3D LV system. In fact, the equations for {X 2 , Y 2 ,Z 2 } are just an integrable 
deformation of the LV equations provided that the constraint X\ = is imposed, and the perturba- 
tion terms contain both the {X 2 , Y 2 ,Z 2 } variables as well as the {Yi,Z\} ones. Complementarily, 
the equations for {Xi ,Y±,Zi} are another integrable deformation the of LV equations under the 
constraint X 2 = 0. Obviously, this dynamical twisting comes indeed from the explicit form of 
the coproduct map AC 2 ) or, in geometric terms, from the fact that the group law that defines the 
coproduct is non-abelian. 

Higher dimensional generalizations of this construction are straightforwardly obtained from the 
definition of the iV-th coproduct map A^ N \ By following the same procedure as in the N = 2 case 
-see (I2.13P and (12. 15ft -. the iV-th coproduct map is obtained by computing the product of N group 
matrices Mn • Mj\r-i M 2 • Mi and reads 



N 



N 



N 



A w (x) =nx i a w (y) = 5> n x f a(n) w = e z > n x 



N 



N 



-7/0 



(3.5) 



i=l 



i=l j=i+l 



2=1 j=i+l 



which is again, by construction, a Poisson homomorphism for the iV-th tensor product of the LV 
bracket V, as it can be explicitly checked. Now, if we define the Hamiltonian T-L^ as the iV-th 
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coproduct of H 



H {N) = 01 A ^ (X) + a 2 (F) + a 3 (Z) + 61 A^' (log X) + b 2 A^ (log Y) + 63 A W (log Z) 

(3-6) 

we get a set of 3N Hamilton equations for the variables {X%, Y^, Z^}, that we do not write explicitly 
here for the sake of brevity. As expected, this system generalizes the N = 2 pattern previuosly 
described: for each fixed subset of three variables {X}~, Y^, Z^} we recover an integrable deformation 
of the LV system provided that the (N — 1) constraints Xj = for all j ^ k are imposed. These 
constraints reflect again the 'twisting' between the ./V different sets of perturbed LV equations that 
arise as a consequence of the Poisson invariance of the system under the iV-th coproduct map. We 
recall that the generic properties of the dynamical embedding of a (N — l)-th coproduct system 
into the corresponding iV-th coproduct one were discussed in [24J in terms of the so-called 'cluster 
variables'. 

Finally, by working out in detail the general formalism of Poisson coalgebra integrability given 
in [20], it can be proven that the set of integrals of the motion in involution for this 3 iV-dimensional 
system are the Hamiltonian together with the N + (N — 1) functions given by the Casimirs on 
each group element and the r-th coproducts (r = 2, . . . , N) of the abstract Casimir C. Explicitly, 





—a 




C (r) ._ A M(£) 


= [A^(X) 




r 












.4 = 1 





i = 1 JV 



e« n 

i=l j=i+l 



e* n x, 

i=l j=i+l 



■<y/a 



(3.7) 



r = 2,...,N (3i 



where we have used the superscript (r) to denote the comultiplication map between M and the 
tensor product of r copies of M. 



4 Integrable perturbations from a deformed Poisson-Lie group 

This group-theoretical interpretation of the LV equations (|2.4p raises also the natural question 
concerning the existence of deformations of the Poisson structure V (12. If) that could be also invariant 
under the coproduct map A^ 2 ). Equivalently, this is just the problem of finding other PL structures 
on the group G, whose existence would provide new integrable deformations of the LV equations. 

This question can be answered in the affirmative, and the following result can be obtained 
through direct computation: 

2. 7 

Proposition. The most generic quadratic Poisson structure in {X, X? , X~a ( Y, Z, 1} and for 
which the comultiplication A^ 2 ) is a Poisson map is given by the brackets 

{X,Y} = aXY + 6X(l-X0) 

{X,Z} = (3XZ + eX(l - X~i) (4.1) 

{Y,Z} = jYZ + j-Y + ^Z + ^—h-XP »). 

p a pa V / 

We shall call this Poisson bivector as Vs, e - Moreover, the Casimir function for Vs, e is found to be 



8(l-Xf>) + ctY a e(X% -1) + /3ZXZ 



(4.2) 



Consequently, we can say that ('Pa^A^ 2 )) is a multiparametric PL structure on G whose limit 
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5, e — > leads to the LV Poisson structure V. 

With this result in mind, an integrable (5, e)-deformation of the LV equations is straightforwardly 
obtained by considering again the same Hamiltonian (|2,3p and the new Poisson bracket Vs, e - In 
this way we obtain the perturbed LV equations 



X = X [aa 2 Y + f3a 3 Z + (ab 2 + f3b 3 )} 

+5X (l - xV} (a 2 + y J + eX (l - X~l 



&3 

a3+ z 



Y 



Y [-aaiX + ja 3 Z + (-ah + 763)] +5 {X* - l)(aiX + h) + 1 (a 3 Z + b 3 ) 

a 



' P 



1 - XTi~* 



(4.3) 



Z [-P<nX - ja 2 Y + (-762 - Ph)] 

1 



Z 



a 2 + 



Y 



+e 



[x- 



l^aiX + h 



(o 2 Y + b 2 



+ 



Se'y 
aP 



(XI 



1) aa + 



Y 



whose integrability is preserved despite the presence of both polynomial and rational perturbation 
terms, since both H and the deformed Casimir C$, e are integrals of the motion in involution for the 
system. 

Obviously, the (5, e)-integrable deformation of the iV sets of twisted LV equations given in 
Section 3 can be obtained by using Vs,e (instead of V) and by making use of the same iV-th coproduct 
map (|3.5p and Hamiltonian (|3.6p since Vg e is invariant under the same group multiplication. Again 
we do not present the explicit form of the set of 3N equations here, but we would like to stress that 
in this deformed case the following (2N-1) integrals of the motion for the system can be explicitly 
obtained from the coalgebra integrability approach: 



S[l-Xf 



+ aYi 



AW(C) = \S (l - A^(X)i) + aA^(Y) 



e AT* -1 +pZiX, 



1,2, 



./3 r 



(A) « - 1) + /3AW(Z)A^(A)« 



Indeed, these are the appropriate deformations of the coalgebra integrals (|3.7p and (13. 8p , 



(4.4) 

r = 2,..., 
(4.5) 



Concluding remarks 

Summarizing, we have presented a new Poisson-Lie interpretation of the quadratic Poisson structure 
that underlies the integrability of the 3D LV system and we have shown that this additional group 
invariance property is useful for getting new integrable deformations of the LV system through 
two different procedures. The first one makes use of the invariance of the Poisson bracket under 
the product of N group elements, and gives rise to a nested set of N deformed LV systems. The 
second one provides a two-parametric familly of integrable deformations of the 3D LV equations by 
finding the most generic deformation of the quadratic Poisson bracket that is compatible with the 
group multiplication. Therefore, this quantum-group algebraic machinery (we recall that Poisson- 
Lie groups are just the 'semiclassical' limit of quantum groups [H]) turns out to be again useful in 
a classical integrability context. 

Concerning applications, it would be certainly intreresting to find specific LV models in which 
the (5, e)-perturbation terms appearing (14. 3D could be dynamically meaningful, and in this context 
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the presence of rational perturbation terms could give rise to dynamical regimes with saturation 
rates. Also, the generalization of the Poisson-Lie group approach here presented to the case of 4D 
(and higher dimensional) LV systems is under investigation. 

Acknowledgements 

This work was partially supported by the Spanish MICINN under grant MTM2010-18556 and by 
INFN-MICINN. 

References 

[1] A.J. Lotka, Elements of Mathematical Biology, (Dover: New York) (1956). 

[2] V. Volterra, Lecons sur la Theorie Mathematique de la Lutte pour la Vie, (Gauthier Villars: 
Paris) (1931). 

[3] M. Plank, J. Math. Phys. 36, 3520 (1995). 

[4] B. Hernandez-Bermejo, V. Fairen, Mathematical Biosciences 140, 1 (1997). 
[5] T.C. Bountis, M. Bier, J. Hijmans, Phys. Lett. A 97, 11 (1983). 
[6] L. Cairo, M.R. Feix, J. Goedert, Phys. Lett. A 140, 421 (1989). 

[7] B. Grammaticos, J. Moulin-Ollagnier, A. Ramani, J.M. Strelcyn, S. Wojciechowski, Physica 
A 163, 683 (1990). 

[8] L. Cairo, M.R. Feix, J. Goedert, J. Math. Phys. 33, 2440 (1992). 

[9] L. Cairo, J. Llibre, J. Phys. A: Math. Gen. 33, 2395 (2000). 
[10] K. Constandinides, P.A. Damianou, Reg. Chaotic Dyn. 16, 311 (2011). 
[11] Y.T. Christodoulides, P.A. Damianou, J. Nonlin. Math. Phys. 16, 339 (2009). 
[12] Y. Nutku, Phys. Lett. A 145, 27 (1990). 
[13] P.A. Damianou, Phys. Lett. A 155, 126 (1991). 
[14] H. Giimral, Y. Nutku, J. Math. Phys. 34, 5691 (1993). 
[15] M. Plank, Nonlinearity 9, 887 (1996). 
[16] E.H. Kerner, J. Math. Phys. 38, 1218 (1997). 
[17] B. Hernandez-Bermejo, V. Fairen, J. Math. Phys. 39 6162 (1998). 
[18] Puyun Gao, Phys. Lett. A 273, 85 (2000). 

[19] V. Chari, A. Pressley, A Guide to Quantum Groups, (Cambridge University Press: Cambridge) 
(1994). 

[20] A. Ballesteros, O. Ragnisco, J. Phys. A: Math. Gen. 31, 3791 (1998). 

[21] A. Ballesteros, A. Blasco, F.J. Herranz, F. Musso, O. Ragnisco, J. Phys: Conf. Ser. 175, 
012004 (2009). 

9 



[22] F. Musso, J. Phys. A: Math. Theor. 43, 434026 (2010). 
[23] F. Musso, J. Phys. A: Math. Theor. 43, 455207 (2010). 
[24] A. Ballesteros, O. Ragnisco, J. Phys. A: Math. Gen. 36, 10505 (2003). 



10 



